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O ■ Abstract 

We present preliminary results of searches for exclusive B^ decays to tt^tt^tt^ among 22.7 million bb 
^ . pairs collected by the BABAR experiment from electron-positron collisions near the T(4S) resonance. 

We measure B{B^ — > /o^vr ) = (28.9 it 5.4 ± 4.3) x 10~^, and find no evidence for the presence of 
any other decay mode in the 'k~^'k~'k^ Dalitz plot. Upper limits are determined for the branching 
,^ . fractions of B^ — > p'^vr'^, non-resonant B^ decays to vr^vr^vr'^ and of several discrete regions of 

■K^TT^TT^ phase-space. We also measure the direct CP-violating asymmetry between the rates of 
untagged p^Tx~ and p^vr^, finding no significant evidence for an effect. 
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1 Introduction 

We describe the measurement of branching fractions for B^ decays]^ to the final state ■K~^'ir~Tr^. We 
include here those decays which proceed through a resonant, quasi-two body intermediate state, as 
well as non-resonant three-body decays with which they interfere. 

The intrinsic interest in these modes stems from their potential use for measuring the angle a of 
the unitarity triangle |||]. Such measurements would exploit interference between the B^ -^ p^vr 
modes and the colour-suppressed B^ — > p^rr^ mode. However, the former modes were discovered 
only in 1999 Q, and the latter mode remains undiscovered. Furthermore, the proposed measure- 
ment of a assumes that only the lowest-mass p mesons contribute to the tt^-k"!!^ Dalitz plot, an 
assumption that needs to be tested. It has also been suggested HI, §] that resonances with quantum 
numbers other than those of the p might make significant contributions to the rate in the Dalitz 
plot, and it is important also to search for these. In the analyses presented here, we improve on 
earlier measurements of B^ — > p^vr and search for B^ — > p^ir^ modes containing other resonances, 
and non-resonant B^ -^ tt^tt^tt^ decays. 

2 The BaBAR Detector and Dataset 

The data used in these analyses were collected with the BABAR detector at the PEP-II storage 
ring. The BABAR detector, described in detail elsewhere Q, consists of five active sub-detectors. 
Surrounding the beam-pipe is a silicon vertex tracker (SVT) to track particles of momentum less 
than ~120 MeV/c and to provide precision measurements of the positions of charged particles of all 
momenta as they leave the interaction point. A beam-support tube surrounds the SVT. Outside this 
is a 40-layer drift chamber (DCH), filled with an 80:20 helium-isobutane gas mixture to minimize 
multiple scattering, providing measurements of track momenta in a 1.5 T magnetic field. It also 
provides energy-loss measurements that contribute to charged particle identification. Surrounding 
the outer circumference of the drift chamber is a novel detector of internally reflected Cherenkov 
radiation (DIRC) that provides charged hadron identification in the barrel region. This consists of 
quartz bars of refractive index ~1.42 in which Cherenkov light is produced by relativistic charged 
particles. This is internally reflected and collected by an array of photomultiplier tubes, which 
enable Cherenkov rings to be reconstructed and associated with the charged tracks in the DCH, 
providing a measurement of particle velocities. Outside the DIRC is a CsI(Tl) electromagnetic 
calorimeter (EMC) which is used to detect photons and neutral hadrons, and to provide electron 
identification. The EMC is surrounded by a superconducting coil which provides the magnetic 
field for tracking. Outside the coil, the flux return is instrumented with resistive plate chambers 
interspersed with iron. 

The data sample used for the analyses contains 22.74 million BB pairs Q, corresponding to 
20.7 fb~ taken on the T(4S) resonance. In addition, 2.6 fb~ of data taken off-resonance have been 
used to validate the contribution to backgrounds resulting from e~^e~ annihilation into light qq pairs. 
These data have all been processed with reconstruction software to determine the three-momenta 
and positions of charged tracks and the energies and positions of photons. Refined information on 
particle type from the various sub-detectors described above is also provided, and is used in particle 
identification algorithms in the analyses. 



^Charge conjugate initial and final states are assumed everywhere unless otherwise stated. 







3 Analysis Method 

3.1 Dalitz Plot Decomposition for 7r"'"7r~7r 

Data for the B^ -^ tt^tt^tt^ final state can be represented on a Dalitz plot. A complete picture 
of the Dalitz plot structure can be obtained only by performing an amplitude analysis, but the 
simplest such analysis uses eight parameters |^], and would be difficult with the relatively low 
statistics currently available. Previous analyses of such decays in this and other experiments 0, 1^, H 
have focused predominantly on resonant quasi two-body decays to /9(770)7r, in which the data are 
localised in bands within the Dalitz plot. We follow this approach, except that our current statistics 
are sufficient to allow also searches for higher resonances populating different bands in the Dalitz 
plot. We therefore sub-divide the Dalitz plot into several distinct regions, each of which is chosen 
to be sensitive primarily to a single resonance, such as the p{770), /9(1450) and /o(400— 1200) (also 
called a), and we perform an independent search in each. In addition, we consider the middle part 
of the Dalitz plot, which is sensitive to non-resonant decays. We have found that the measured 
background level varies dramatically across the Dalitz plot and our approach allows the analysis 
for each resonance to be optimised independently, taking into account the local background level. 
The seven regions are indicated in Fig. |l[ The regions are defined using selection criteria for the 
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Figure 1: The branching fraction into the Dalitz plot is measured in separate regions which are sen- 

III: 5° -^ p'^TT^, IV: B^ -^ p'^tt'^, V: 
TT^TT^TT^ at high mass. 
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sitive to the following modes: I: B^ — > p^vr^, II: B 

B^ — > charged scalar and vr^, VI: B^ -^ neutral scalar and n^, VII: B^ 

The dashed line along the diagonal represents the CP-symmetry axis. 



invariant mass of vrvr-pair combinations, m.,^Tj, and the pair helicity angle. Oh, defined as the angle 
between the direction of one of the pions and the direction of the parent B meson candidate, 
computed in the vrvr-pair rest frame. The region definitions are summarised in Table [l|. 



Table 1: Regions of the B^ — > tt'^tt tt^ Dalitz plot: the regions are kinematically-defined by the 
and I cos {9h)\ selection criteria. Any possibly dominant resonance is noted for each region. 



m 



Region 


Putative dominant 
resonance 


m,r7r(GeV/c^) and cos (0//) selection 


I 


p±(770) 


0.52 < m^+^o < 1.02 and | cos {9+")\ > 0.25 
0.52 < m^-^o < 1.02 and cos (e^f)] > 0.25 


II 


^^^(770) 


0.57 < m^+^- < 1.02 and | cos (6*^)1 > 0.3 
and m^+^o > 1.2 and m^-^o > 1.2 


III 


/9±(1450) 


m^+^o and m^-^o and m^+^- > 1.2 

and ((m^+^o < 1.8 and | cos (6'+°)| > 0.25) 

or (m^-^o < 1.8 and | cos {9^^)\ > 0.25)) 


IV 


p"(1450) 


1.2 < m^+^- < 1.8 and cos 6*;^" > 0.3 
and m^+^o and m^-^o > 1.9 


V 


charged scalar 


(0.52 < m^+^o < 1.8 and | cos {ep')\ < 0.25) 
or (0.52 < m^-^o < 1.8 and | cos {e]f)\ < 0.25) 


VI 


/'^(400 - 1200) 


0.57 < m^+^- < 1.8 and | cos (6lJ~)| < 0.30 


VII 


high mass(> 1900) 


m'^+^o and rn-^_^o and m^_|_^_ > 0.13 m'^ 



3.2 Candidate Selection 

Charged tracks are required to have a momentum less than 10 GeV/c and a transverse momentum 
greater than O.lGeV/c. They are required to have at least 20 hits in the DCH and to originate 
close to the beam-spot. In addition, vr candidates are selected as those charged tracks which 
are unlikely to be kaons, based on Cherenkov angle information from the DIRC combined with 
energy-loss information from the DCH. 

Photon candidates are identified in the calorimeter as deposits of energy unassociated with 
charged tracks, n^ s are reconstructed by combining pairs of photon candidates and requiring 
that the invariant mass of the resultant candidate is between lOOMeV/c^ and 160MeV/c^. Photon 
candidates used in vr'' reconstruction are required to have a minimum energy of 30MeV and a 
shower shape consistent with a photon. 

Reconstruction of B candidates is accomplished by forming all combinations of 7r+7r~7r'^ can- 
didates in each event and applying a loose quality requirement of x^ < 200 to the fitted vertex 
for the 7r"'"7r~ pair. The B candidates are required to satisfy kinematic constraints appropriate 
for B mesons. We use two kinematic variables [^ for this: iriES = v (2'^ + Po " Pb)^/^o ~ ?*£' 
where the subscripts and B refer to the e^e~ system and the B candidate, respectively; and 
AE = E^ — \/s/2, where E'^ is the B candidate energy in the center-of-mass frame. For signal 
events, the former has a value close to the B meson mass and the latter should be close to zero. 

The selection criteria in each region are optimized for sensitivity to the branching fraction. The 
criteria that are allowed to vary are the mass requirements on the vr'^ and p{770) candidates, cos 6h, 
and the requirements on event shape discussed in the following section. 
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3.3 Background Suppression and Characterisation 

Charmless hadronic modes suffer very large amounts of background from random combinations 
of tracks, mostly from light quark and charm continuum production. Such backgrounds may be 
reduced by selection requirements on the event topologies computed in the T(4S) rest frame. We use 
cos 9t, the cosine of the angle 9t between the thrust axis of the B meson decay and the thrust axis 
of the rest of the event. For continuum-related backgrounds, these two directions tend to be aligned 
because the reconstructed B candidate daughters generally lie in the same jets as those in the rest 
of the event. By contrast, in B events, the decay products from one B meson are independent of 
those in the other, making the distribution of this angle isotropic. In consequence, requiring that 
this opening angle be significant provides a strong suppression of continuum backgrounds. 

Other event shape variables also help to separate signal and background. We form a linear 
combination of 11 variables in a Fisher discriminant J- P]. The coefficients for each variable are 
chosen to maximize the separation between training samples of signal and background events. The 
variables contained in T are: 

• the cosine of the angle between the B momentum and the beam axis; 

• the cosine of the angle between the thrust axis of the B candidate and the beam axis; and 

• the summed momentum, in nine cones coaxial with the thrust axis of the B candidate, of all 
detected particles in the event that are not associated with the B decay. 

Despite the power of such topological variables to reduce the combinatorial backgrounds, most 
of the modes we have searched for continue to suffer significant background levels. Even after 
stringent selection criteria have been applied, it is necessary to make a background subtraction 
to isolate a signal or set an upper limit. In order to do this, the background in the signal region 
is estimated from the number of events in a sideband region, located near the signal region in 
the m-p^Q-^E plane, and extrapolated into the signal region. The shape of the distribution of the 



background as a function of rriES is parameterised according to the ARGUS function |10|, and 
is measured using on-resonance data slightly displaced from the signal region in the /S.E variable 
(0.1 < \/S.E I < 0.25). We define IZ to be the ratio of the number of candidates in the signal region 
to the number in the sideband region. Off-resonance data yields a consistent ratio and is averaged 
with the result measured on-resonance. We measure the value of TZ independently for each region 
of the Dalitz plot, and we find that it is significantly dependent on it. 

3.4 Branching Fraction Analysis 

The branching fractions are calculated according to 

where A^i is the number of candidates in the signal region (SR) for on-resonance data; A'^2 is the 
number of candidates in on-resonance data observed in the sideband region (GSB), so that IZN2 
is the estimated number of background candidates in the signal region; N^^ is the number of BB 
pairs produced and e is the signal efficiency. The GSB is specified by 5.21 < rriES < 5.25GeV/c^, 
\IS.E — {AE) I < 0.1 GeV, where {AE) is the mean value of AE for B meson decays found in 
calibration samples. 
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The numbers and distributions of candidates within the signal region remained unknown to us 
until all aspects of the analysis were finalized. The final selection criteria were chosen to maximize 
the sensitivity to the signal. Once chosen, neither the background description nor the cuts were 
changed. This procedure was followed independently for each channel. 

For the signal efficiency in Eq. Q, we used simulated signal events and the same selection 
criteria as used for the data. The efficiencies due to vr'^ reconstruction, particle identification, and 
AE selection criteria are determined from independent control samples derived from the data. 

We have adopted two slightly differing approaches to the definition of selection efficiencies in 
Eq. m. For the B^ — > p^ir^ and B^ — > p^i:^ measurements, we simulate the p resonances according 
to a non-relativistic Breit-Wigner distribution within the phase-space and with a cos^ 9h helicity 
distribution. For the non-resonant B^ -^ tt^tt^tt^ measurement, we simply assume a uniform 
distribution in phase-space, and ignore interferences. In these cases, we take as the denominator in 
the definition of the selection efficiency, the number of events generated in the full Dalitz plot, so that 
the loss of events outside our analysis cuts is compensated in Eq. 1^. This provides a conventional 
branching fraction, having assumed a signal distribution, and ignoring interferences. By contrast, 
for the measurements in the remainder of the Dalitz plot, we have no a priori model of the dynamical 
distribution of events in phase space or of the feedthrough between different modes (other than of 
B^ — > p^vr^, whose rate is known), so we simulate the events according to a uniform distribution 
in phase space. We note that this introduces no systematic error, provided that the efficiency is 
flat over the region. We have verified that this is the case within the statistics available to us (the 
efficiency away from the kinematic boundaries, projected onto either Dalitz plot axis is consistent 
with a linear dependence, and the slope has been shown to be zero within errors of magnitude 
±3 X 10^^ GeV^ ). In these cases, we take as our denominator for the efficiency calculation, only 
those events generated within the boundary of the defined region. In this way, events lying outside 
this boundary are not considered part of the signal. For these latter "topological" branching 
fractions, the measurement should be considered the total rate within the defined boundaries, not 
just that associated with any single resonance. If a signal were found, more work would be needed 
to determine the detailed dynamics. 

3.5 Charge Asymmetry in B^ — > p^ir 

In principle, there are four decay rates that we could measure for the generic B^ — > p^vr^ mode: 
we define 

r(50 ^ P+7T-) = Tp^ ; r(50 ^ p-7T+) = r^p (2) 

and their CP conjugates 

r(_BO ^ p-TT+) = Tp^ ; r{B^ ^ p+TT-) = T^p (3) 

respectively. In the general, CP-violating case, all four are different, and there are two observable 
direct CP violations: 

and A^p = (F^p - r^p) / 0. (4) 

In the present analysis, we have not implemented tagging, but we have divided our B^ -^ p^vr 
sample into two sub-samples, containing respectively the events with p+vr" and those with p'-K^ . 
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Using these, we may form the charge asymmetry: 



•^ " N{p+7r-) + N{p-TT+) ~ •^P'^^^ + ■^'^'' ^^^ 

where ^det is any detector-induced charge bias and 

V^^ p-TT I J- np) + (J- p7r + i -j^pj 

The numerator of Eq. g is simply the difference of the two direct CP violations in Eq. |^. If it 
is non-zero, it confirms direct CP-violation in at least one of the decays, Eqs. |2| and ^. We may 
assume that systematics due to background subtraction and luminosity cancel in the numerator of 
Eq. ^. Detector charge biases from tracking and from particle identification cuts were measured 
in an independent study with control samples of identified D"^ -^ ^s'^^ events. Single it charge 
asymmetries were limited to < 2%, so we conservatively set a ±2% systematic uncertainty on our 
asymmetry measurement. 

4 Treatment of Experimental Uncertainties 

The systematic errors associated with branching fraction measurements arise from uncertainties in 
the background subtraction, in the overall signal efficiency and in the overall normalisation N^^. 

In general, there are two contributions to the background subtraction: continuum backgrounds 
and -B-related backgrounds. The estimate of the continuum background is given as the product 
of the number of events in the GSB and the factor TZ, which is the estimated ratio of the number 
of background events in the GSB to that in the signal region, and is "measured" with sidebands 
and off-resonance running as described above. The number of events observed in the GSB has 
only a statistical error, which is taken into account mode-by-mode. The factor 7^ has a systematic 
error given by the error in the fitted ARGUS shape parameter, ^. This is taken into account 
independently for each mode. 

Possible sources of -B-related backgrounds include events with low-multiplicity decays to charm, 
which have been excluded explicitly by cuts in the Dalitz plot, and other charmless decays. The 
latter have been dealt with by counting the numbers found in the signal region in a dedicated 
charmless B decay "cocktail" simulation, scaling them according to the relative luminosities in data 
and simulation samples, and subtracting the appropriate numbers of events from the numbers in the 
signal region. We include relevant statistical errors from the simulation sample, and conservatively 
set the systematic error in this subtraction to be equal to the number subtracted. 

We measure the overall signal efficiency e with signal Monte Carlo simulation, and some un- 
certainty arises from limited signal mode statistics. The fractional error from this source is 2.7% 
in the B^ — > p^ir mode. In the other modes, this error is negligible compared to other sources of 
uncertainty, which are generally larger than the observed signals. The accuracy of the simulation is 
subject to systematic uncertainties including the tracking efficiency, calorimetric shower efficiency, 
particle identification efficiency and the interaction of selection requirements with the resolution 
accuracy of the simulation. At present the dominant uncertainties are in the particle identification, 
TT^ reconstruction, and AE and event-shape modelling in the simulation. 

A discrepancy between the measured track-finding efficiency, and that obtained in our simu- 
lation is taken into account with a momentum-dependent weight per track. The uncertainty in 
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the determination of these weights leads to an uncertainty in the tracking efficiency of 1.2% per 
track, added coherently for both tracks used in the analyses. For vr^ reconstruction efficiency, we 
provide an energy-dependent weight, which is averaged over all accepted candidates in our signal 
simulation samples to provide a sample-dependent weight. We conservatively assign a systematic 
uncertainty of 5% per vr*^ to this weight. The vr^ particle identification efficiency was measured in 
data with a vr control sample. It was found to be significantly smaller than that predicted in the 
simulation, and this difference is taken into account with a factor i?™'"'" = 0.948 it 0.018, per vr^, 
added coherently for both pions. 

A control sample of B~^ -^ D p^ ^ K^tt^tt^tt^, which has similar kinematics to our modes, 
has been used to show that the simulation models the AE shape of such signals quite well. However, 
due to statistical limitations in this method we ascribe a systematic uncertainty of 5% to the signal 
efficiency from our simulation, to account for possible resolution differences on the order of 5 MeV 
in AE. Such uncertainties in the rriES selection efficiency are negligible by comparison. For event 
shape cuts, the uncertainty on the efficiency has been taken as 5%. 

The overall normalisation, N^^, comes from a dedicated "S-counting" study, which was made 
to determine the number of B mesons in the data samples. The final systematic error on this was 
determined to be 1.6%. 

The overall systematic uncertainty is the sum in quadrature of the contributions from all sources. 
Table ^ shows a summary of sources and magnitudes of systematic errors. 



Table 2: Summary of Systematic Error Sources and Magnitudes. 



Source of Uncertainty 


Uncertainty 
as % of Result 


B" ^ pT 


^± 


Other modes 


Background subtraction: 

A 

S-related Backgrounds 




5.2 

5.8 


> 10 
>20 


sub-total 


7.8 


> 30, dominant 


MC Statistics 

MC Efficiency Corrections 


2.7 
10.5 


negligible 
negligible 


_B-counting 


1.6 


1.6 


Total 


14.5 


> 30 



5 Results 

Our preliminary measurements of branching fractions are summarized in Table 3. In calculating 



the upper limits, we have used the classical method outlined in [11|. In order to take account of 
the systematic errors, we have reduced the background estimate and the efficiency by one standard 
deviation (systematic) before making the calculation. 

The new value for the B^ — 
ments. Our upper limit for B^ 



p^TT branching fraction is more precise than previous measure- 



p^vr^ is consistent with Standard Model expectations |12]. Our 
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Table 3: Results for the Svr final state analyses. We aggregate the first two columns to obtain 
the branching fraction for B^ -^ p ir^ (the asymmetry is calculated independently.) 



oi 



Quantity 


5° ^ p+TT- 


5° ^ p-7r+ 


B^ ^ p\^ 


B° -^ 7r+ TT- vrO (NR) 


cos 9t 


<0.75 


<0.75 


<0.65 


<0.85 


Fisher 


<0.75 


<0.75 


<0.65 


< 1.1 


TfO mass GeV/c2 


0.110-0.160 


0.110-0.160 


0.115-0.155 


0.115-0.155 


events in SR 










On-res data 


121 


118 


27 


41 


Off-res data 


11 


12 


5 


3 


events in GSB 










On-res data 


590 


540 


142 


362 


Off-res data 


74 


74 


24 


46 


A 


0.128 ±0.005 


0.128 ±0.005 


0.132 ±0.012 


0.116 ±0.006 


qq BG in SR 


75.5 ±3.1 ±3.0 


69.1 ±3.0 ±2.7 


18.7 ±2.2 ±1.7 


42.0 ±3.1 ±2.0 


bb BG in SR 


2.7 ±1.6 ±2.7 


2.7 ±1.6 ±2.7 


2.2 ±1.5 ±2.2 


3.2 ±1.8 ±3.2 


Signal 


42.8 ±11.5 ±4.0 


46.2 ±11.4 ±3.8 


6.1 ±5.8 ±2.8 


-4.2 ±7.3 ±3.8 


Signal Efficiency 


0.135 ±0.016 


0.074 ± 0.009 


0.075 ± 0.010 


Stat, signif. (a) 


5.01 


0.96 


N/A 


B (xlO-6) 


28.9 ±5.4 ±4.3 


3.6 ±3.5 ±1.7 


N/A 


90% CL limit (xlO-6) 


N/A 


< 10.6 


<7.3 



upper limit of 7.3 x 10~^ for non-resonant Tr+vr^vr^ decay ignores interference effects, an approacli 
similar to that adopted by CLEO ||l^ for non-resonant B~^ decays to three charged pions. 

In addition, we have made a preliminary measurment of the CP asymmetry defined in Eq. (P) 
and find 

Aphys = -0.04 ± 0.18 ± 0.02, (7) 

consistent with zero. This does not however indicate that there is no direct CP violation in 
B^ -^ p^TT^ decays, as the measured quantity is the difference between two direct CP violations, 
and we cannot exclude the possibility of a (coincidental) cancellation of the two. 

The signal for B^ -^ p^ir^ can be seen in the distributions of niEs and ^E shown in Fig. |2|. 
The results for the topological branching fractions are given in Table 4. The existence of a slight, 




Figure 2: mss and /S.E distributions for B^ -^ p tt^. The signal region requirement was made on 
the orthogonal variable in each case. 



though not statistically significant excess in Region III could be a hint of structure due to higher 
resonances, such as the p'^(1450). This would have implications for the measurement of a [p!^ . 
Further data are necessary to confirm whether this effect is real. 

Region VI is sensitive to the scalar resonance /o(400 — 1200), proposed as a significant contri- 
bution in this region Q. Our sensitivity has not yet reached that required to test this possibility. 

We can get an impression of the whole Dalitz plot from Fig. y, which shows the full Dalitz 
plot area for various samples. Fig. |3| a) shows the distribution of the data in phase space after the 
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Table 4: Results for the Svr final state analyses. 



-a 



Quantity 


Region III 


Region IV 


Region V 


Region VI 


cos^T 


<0.75 


<0.65 


<0.75 


<0.55 


Fisher 


<0.75 


<0.70 


<0.75 


<0.40 


7r° massGeV/c^ 


0.110-0.160 


0.115-0.155 


0.110-0.160 


0.115-0.155 


events in SR 










On-res data 


75 


8 


44 


6 


Off-res data 


4 


2 








events in GSB 










On-res data 


390 


80 


268 


22 


Off-res data 


54 


6 


38 


2 


A 


0.131 ±0.006 


0.134 ±0.011 


0.132 ±0.009 


0.142 ±0.017 


qq BG in SR 


51.1 ±3.7 ±2.4 


10.7 ±1.7 ±0.8 


35.4 ±3.1 ±2.3 


3.1 ±0.9 ±0.4 


bb BG in SR 


6.5 ±2.5 ±6.5 


2.0 ±1.4 ±2.0 





3.2 ±1.8 ±3.2 


Signal 


17.4 ±9.7 ±6.9 


-4.7 ±3.6 ±2.2 


8.6 ±7.3 ±2.3 


-0.3 ±3.2 ±3.2 


Signal Efficiency 


0.150 ±0.021 


0.087 ±0.018 


0.150 ±0.023 


0.067 ±0.014 


Stat, signif. (a) 


1.83 


N/A 


0.40 


N/A 


B (xl0~6) 


5.1 ±2.9 ±2.2 


N/A 


2.5 ±2.1 ±0.8 


N/A 


90% CL limit (xlO-6) 


< 11.3 


<2.7 


<6.1 


<5.2 




% 30 
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Figure 3: Dalitz plots for B^ — > tt^tt tt^ after i?*^ 



p^vr^ 



analysis requirements (except m^Tr and 



\cos9h\)'- (a) data SR, (b) data GSB, (c) Monte Carlo SR, (d) background-subtracted data SR. 
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analysis cuts are applied, and Fig. H b) shows the background distribution measured in the GSB. 
Figure ^ c) shows an example Monte Carlo signal and Fig. ^ d) shows the background-subtracted 
signal. The helicity structure is clearly visible along the two p -bands in the background-subtracted 
figure. The peak near the origin {ie. low ir^ energy), even after background subtraction, indicates an 
efficiency enhancement in real signal events (also present in the simulated signal sample), in which 
events missing one low-energy photon due to acceptance effects are recovered by combinatorial 
substitution of a different photon. This effect significantly enhances the population of signal events 
at low vr'' energies (only one candidate is accepted per event). 

6 Conclusion 

We have made a number of preliminary measurements of branching fractions into regions of the 
B^ — > 7r"''7r~7r'^ Dalitz plot. In particular, a new measurement of the B^ -^ p^n branching fraction, 
upper limits for the B^ — > p^n^ mode and non-resonant B^ -^ 7r"'"7r~7r'^, and four upper limit 
measurements of topological branching fractions into other regions of the Dalitz plot have been 
reported. We have also made the a preliminary measurement of the direct C-P- violating asymmetry 
between the rates of untagged p"^7r~ and p^vr^, finding no significant evidence for an effect. 
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